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The interplay between spontaneous symmetry breaking in many-body systems, the wavelike na¬ 
ture of quantum particles and lattice effects produces an extraordinary behavior of the chiral current 
of bosonic particles in the presence of a uniform magnetic flux deflned on a two-leg ladder. While 
non-interacting as well as strongly interacting particles, stirred by the magnetic field, circulate along 
the system’s boundary in the counterclockwise direction in the ground state, interactions stabilize 
vortex lattices. These states break translational symmetry, which can lead to a reversal of the 
circulation direction. Our predictions could readily be accessed in quantum gas experiments with 
existing setups or in arrays of Josephson junctions. 

PACS numbers: 03.75.Lm, 05.30.Jp, 37.10.Jk 


A charged quantum particle, due to its wavelike na¬ 
ture [ij, picks up an increment to its phase proportional 
to the magnetic flux piercing the area enclosed by the 
particle’s (closed) path. If the particle hops around a 
plaquette, the accumulated phase is proportional to the 
magnetic flux $ = Ba^, a being the lattice constant. 
Since the phase of the wave-function is defined modulo 
27r, the action of a magnetic field on quantum particles 
in a lattice is periodic: If we define a dimensionless flux 
(j) = 27r$/<I>o with the magnetic flux quantum $o = h/q, 
where h is Planck’s constant and q is the charge, any 
physical quantity A obeys A{4>) = A^cf) + 27r). The mini¬ 
mal extensive lattice on which this behavior can emerge 
is a two-leg ladder (see Fig. [1]). 

In this work, we discuss the intriguing effect of a re¬ 
versal of the circulation direction of the chiral current of 
interacting bosons on a two-leg ladder, due to the spon¬ 
taneous formation of a large unit cell without changing 
the external magnetic field. The key ingredients to real¬ 
ize this effect are, first, the wavelike nature of quantum 
particles defined on a lattice, and second, many-body ef¬ 
fects. 

The basic idea is sketched in Fig. [H For a single pla¬ 
quette, the magnetic fields corresponding to values of the 
flux 0 < (f> < TT produce a ground-state net current with 
a counterclockwise chirality. When one assembles these 
plaquettes into a minimal extensive lattice such as the 
two-leg ladder, naively, one would expect that the local 
currents on individual links along the boundary of the 
ladder add up to produce a net chiral current jc, also 
circulating in counterclockwise direction for 0 < </> < tt. 
Indeed, this is the case for non-interacting particles 0. 
In general, since the chiral current is a ground-state prop¬ 
erty, the relevant unit cell does not need to be the same 
as the unit cell of the underlying lattice. 
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FIG. 1. (Color online) A spontaneous doubling of the unit 
cell leads to an increase of the effective flux as well: (p —>■ 2(j>. 
Under the conditions described in the text, the chiral current 
jc (indicated by the arrows) reverses its direction. 


If the unit cell is doubled, then this will result in a 
doubling of the effective flux = 2(j), piercing the 
enlarged unit cell of the ground state. The chiral cur¬ 
rent jc is an odd function of flux and periodic, thus 
jc{4>) = jc{4> — 27r) = —jc(27r — <f>). As a consequence, 
for values of the flux 7r/2 < cp < tt defined with re¬ 
spect to the original cell, the effective flux is in the do¬ 
main £ (~^j0) modulo 27r, which one would also 
obtain if the orientation of the magnetic field was in¬ 
verted. Here, we show that such an increase of the unit 
cell arises as a consequence of the spontaneous breaking 
of discrete lattice-translation symmetry in vortex-lattice 
(VL) phases. 

We study the single-band Bose-Hubbard model defined 
on a two-leg ladder lattice with flux cp per plaquette [see 
Fig.EKa)]: 
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aj ^ creates a boson on the lower (£ = 1) or the upper site 
[(. = 2) of the r-th rung, n^^j. = aj and L is the num¬ 
ber of rungs. The hopping matrix elements between the 
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FIG. 2. (Color online) (a) Sketch: Two-leg ladder with 
L rungs and a uniform flux (j) per plaquette. (b)-(d) Local 
currents and onsite densities in (b) the Meissner phase, (c) 
the VLi/3-SF and (d) the VL 1 / 2 -SF [U = 2J, = 1.6J, 

p = 0.8 and (b) 0 = O.Ott, L = 10, (c) (ji = O.Stt, L = 120 and 
(d) (j) = O.Ott, L = 120]. The length of the arrows encodes 
the strength of local currents. The size of the circles and the 
intensity of the shading are proportional to the onsite density. 
In (c), a Meissner-like region is indicated by the dashed line. 
The chiral current is reversed in (d), but not in (b) and (c). 
The value of Jj_/ J = 1.6 was chosen since two VLs are realized 
for these parameters. 


nearest-neighbor sites along the ladder’s legs and rungs 
are denoted by J and Jj_, respectively, and U is the repul¬ 
sive onsite interaction. The filling is p = N/{2L), where 
N is the number of bosons. We carried out density ma¬ 
trix renormalization group (DMRG) simulations [3|,1^ for 
p <2 and U > J. For the complementary regime of large 
densities p ;§> 1 and weak interactions U Jp, we use a 
mapping to a frustrated XY model and apply a transfer- 
matrix approach (see the supplemental material Q 
for details on both techniques). 

Physical realizations are, for instance, either arrays of 
Josephson junctions of superconducting islands in a mag¬ 
netic field [sMlOll , where the bosonic particles are Cooper 
pairs of electrons with g = 2e, or neutral ultracold bosons 
in optical lattices [ill . 121 in the presence of artificial 
gauge fields, subject to s-wave interactions 13 In 
the latter type of experiments, using either superlattices 
[l3| or a synthetic lattice dimension [H-lIll to realize 
ladders, chiral currents have recently been measured. 

For the two-leg bosonic ladder, the existence of 
Meissner-like and vortex phases has been firmly es¬ 
tablished both in the weakly- and strongly-interacting 
regime 1^ 24|. In Fig. [2Kb), obtained from a DMRG sim¬ 
ulation, we depict the typical behavior of local currents 
and particle densities in the low-field Meissner phase. 
The currents are non-zero along the boundary of the two- 
leg ladder, consisting of the edge rungs and the legs of 
the ladder, and vanish quickly on the inner rungs away 
from the first and last rung, producing a net chiral cur¬ 


rent jc = (ErO'i'r - il.r) + {jr=l -jr=L))/^ (§66 ScC. SI 
in [71 for definitions and properties of jc)- The coherence 
of the relative phase of bosons on the two legs induces 
the Meissner phase in bosonic ladders d, [l^ and hence 
a chiral current can also emerge in the Mott-insulating 
regime for fillings with an integer number of particles per 
rung US [ 23 . 

Upon increasing the flux beyond a critical value, the 
system enters into a vortex phase where local currents on 
the inner rungs develop, resulting in a current configura¬ 
tion that is reminiscent of a Meissner phase disordered 
with vortices (Sin, [^ . These vortices interact repul¬ 
sively with each other [3 , yet for generic vortex densities 
Pv, they are distributed in the system without any peri¬ 
odicity. At certain commensurate vortex densities, VLs 
were predicted to form in the ground state @,[ 13 . DMRG 
results for such VL states are depicted in Figs. |2Kc) and 
(d). In the VL superfluid (SF) at vortex density = 1/3 
[Fig- EIc)], the currents on two neighboring plaquettes 
form a complex [surrounded by a dashed line in Fig.|2Kc)] 
that is a mini copy of the Meissner phase (with ’screening’ 
currents circulating around the complex’ boundary and 
vanishing currents on its inner rung) such that on every 
third plaquette, a vortex resides, around which the cur¬ 
rent circulates in the direction opposite to the behavior in 
the Meissner phase. In Fig. He), the chiral current goes 
in the counterclockwise direction. In the VL at = 1/2 
(VL 1 / 2 -SF), obtained for ^ < tt, a vortex sits on every 
other plaquette [Fig. Hd)] and the direction of the chi¬ 
ral current is reversed. The observation of several stable 
VLs is a main result of this Letter. 

For not too low particle densities and flux values close 



FIG. 3. (Color online) jc versus U/J in the proximity of the 
VLs at = 1/2, for (a) p = 0.5 and (b) p — 0.8, (DMRG 
data, L = 160, (j) = O.Ott, J± = 1.6J). The exact value for 
17 = 0 is ~ 0.08J > 0, independently of p. The small stms 
seen in the VSF phase are finite-size effects (see Sec. S2 in [ 3 ). 
The dashed lines serve as guide-to-the-eye for small values of 
U/J in the VSF phase, where for large system sizes (such as 
L=160), it is hard to converge DMRG with respect to both 
the local number of bosons and the DMRG state space [3, [3- 
The dashed lines connect the U — 0 value to the DMRG 
data point for the smallest value of U/J at which we have 
converged results (see Sec. S6 0). plus we indicated the kink 
expected at the VL 1 / 2 -SF to VSF boundary, based on such 
data as shown in Fig. S7 [ 3 ]. 
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to TT, we generically observe a reversal of the chiral cur¬ 
rent tuned by the interaction strength U/J: jc flows in 
a counterclockwise direction along the ladder’s legs for 
U = 0 and 17 —^ oo, yet for a certain range of interaction 
strengths at low temperatures the circulation direction 
of the current is reversed, as if the majority of particles 
swam against the tide. 

The interaction driven chiral-current reversal is evident 
from the DMRG results presented in Figs.[3lja) and (b). 
There, we plot the chiral current as a function of U/J 
for p = 0.5 and p = 0.8, respectively (see 0 for data for 

p= 1 ). 

We clearly numerically resolve a current reversal, 
which is tied to the presence of the VL at = 1/2. This 
state can be realized both in the superfluid (VL 1 / 2 -SF) 
and the Mott-insulating regime (VL 1 / 2 -MI). The VL 1 / 2 - 
SF phase is neighbored by a vortex superfluid (VSF) for 
small values of U/J. Generally, upon entering into the 
vortex-liquid phases from the VL phases, in the ther¬ 
modynamic limit the absolute value \jc\ of the current 
decreases continuously (while for finite L, small steps 
are seen [23), showing a cusp-like behavior at the phase 
boundary [7|. Other transitions that are crossed in Fig. [3] 
do not leave a fingerprint in the interaction dependence 
of jc for L —>■ 00. For p = 0.5, further increasing U/J 
takes the system from the VL 1 / 2 -SF into a VL 1 / 2 -MI. 
Typically, the sign of the chiral current becomes positive 
again on the large U/J side of the VL 1/2 phases. The 
chiral-current reversal is thus robust against the pres¬ 
ence or absence of a mass gap and a variation of density. 
Moreover, the absolute value of the reversed chiral cur¬ 
rent can exceed the U = 0 value. 

The spontaneous symmetry-breaking leading to the 
M-times enlarged unit cell with i/eff = Mcj) € (—7r,0) 
modulo 27r (M = 2 in Fig. |3l) is a vital ingredient for 
obtaining the current reversal. We argue that yet an¬ 
other requirement is that the dominant contribution to 
the current comes from particles with a large wavelength 
experiencing the effective flux. This is exemplified by 
the behavior of the chiral current in the fully gapped 
VL 1 / 2 -MI [see Fig. EKa)]: Inside this state with a dou¬ 
bled unit cell, the chiral current, as a function of C//J, 
passes through zero and returns back to the original di¬ 
rection of rotation. The localization length of bosons in 
the Mott-insulating phase becomes shorter with increas¬ 
ing U/J, restricting the typical wave-lengths of particles 
and as a consequence, less particles see the enlarged unit 
cell with doubled flux. On these grounds, we expect that 
in a given phase with a spontaneously enlarged unit cell, 
the absolute value of the reversed chiral current attains 
its maximum for the smallest U/J in that phase, consis¬ 
tent with the data presented in Fig. El We emphasize 
that the optimal condition for observing the current re¬ 
versal in two-leg ladders is a doubling of the unit cell. 

Interestingly, one could have deduced the existence of 
the chiral-current reversal from the flux-dependence of 



Jj^/J 


FIG. 4. (Color online) (j) versus J_l/J phase diagram at 
U/J = 2 and p = 0.8. See the text for details. 


the ground-state energy shown in Ref. in the weak- 
coupling regime. From the bosonization analysis of 0 , 
valid for Jj_ <C J, one can also obtain the sign change. 
Yet, neither study actually discussed the reversal effect. 
DMRG results indicate that the chiral-current reversal 
driven by interactions also exists on three-leg ladders [i^ . 

The chiral current can also change its sign in fermionic 
ladders [H 0 as a function of ^ S ( 0 , 7 r). Yet, for 
fermions, the effect exists already in the noninteracting 
case as a result of the band structure while for bosons, the 
interaction-induced spontaneous flux increase is crucial 
for the current reversal. 

To identify the regions in parameter space spanned by 
4>, U, p and J_L, in which the VL phases exist, we present 
a representative ground-state phase diagram in Fig. E] for 
UJ J = 2 and p = 0.8 as a function of flux and J_l (see 
[Tj for a density versus U/J phase diagram). The figure 
illustrates the plethora of phases realized in the model: 
the Meissner superfiuid (MSF), a VSF, a VL 1 / 2 -SF, a 
VL 1 / 3 -SF and a biased ladder phase (BLP). The BLP, 
predicted in mean-field theory | 2 l| , spontaneously breaks 
the symmetry between the two legs by imbalancing the 
density 0. The VL 1 / 2 -SF, which leads to the chiral- 
current reversal, is stable down to at least Jx = J/‘^- 

Since the current reversal is connected with sponta¬ 
neous symmetry breaking in the ground state, a rel¬ 
evant question pertains to the effect of temperature. 
To obtain a quantitative understanding of T > 0, we 
consider the limit of large particle densities p 1 
and weak interactions U <C Jp. Using the transfer- 
matrix approach Q, we obtain the results shown in 
Fig. m where the chiral current [Fig. EKa)] and vor¬ 
tex density [Fig. [51(b)] versus flux are presented for dif¬ 
ferent temperatures. Upon reducing temperature, VLs 
appear at pv = 1/2,1/3,2/5,1/4,1/5,.... These have 
an M = 2, 3, 5,4 5, ...-times enlarged unit cell, respec¬ 
tively Sil llil . 0. At zero temperature, the vortex den¬ 
sity, as a function of flux (f>, exhibits the famous devil’s 
0 0 with a plateau for each corn- 


staircase structure 


mensurate value of py. The reversal of the current in 
the VL 3 / 2 “SF phase clearly survives a finite tempera¬ 
ture, which also applies (though for lower values of tern- 
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FIG. 5. (Color online) (a) Chiral current jc/J and (b) vortex 
density versus flux for Jx = J/2 for p ^ 1, ?7 <C Jp. This 
value of Jx is chosen since here, not only the VL at p„ = 1/2 
exhibits a chiral-current reversal at low T, but also those at 
p« = 1/3 and 2/5 (see the insets (c) and (d)). The temper¬ 
ature is ksT — 0.006J (continuous lines) and ksT = O.IJ 
(dotted lines). The dashed lines are for [/ = 0 at T = 0. (e) 

Largest negative value /)/ = min^jc{4‘) of the reversed cur¬ 
rent for Jx = 0.5J (dashed line) and J± = 1.2J (continuous 
line) as a function of ksT/J. See [ 3 ], Fig. S6, for more data. 


peratures) to the current reversal in the vicinity of the 
VL 1 / 3 -SF state for (f) < 27r/3. Thus, the reversal is more 
stable against thermal fluctuations than the underlying 
VLs. The optimal parameters for observing the current 
reversal in the weak-coupling regime are (p ~ O.Qtt and 
Jx — 1.2J where the reversal occurs for ksT < J/2 0, 
as is shown in Fig. [SKe) as well as in Sec. S5 in 0. 

At extremely low temperatures and for the parameters 
of Fig. [5l a reversal of the chiral current is also visible 
in the p„ = 2/5 VL state with M = 5 for ^ ~ 47 r /5 [see 
Fig. EKd)]. Other VLs [corresponding to the regions in 
which plateaus are formed in the p„ = p«((/)) curve in 
Fig. ini[b)] do not induce a reversal of the chiral current 
even at T = 0 since they are stable at values of fluxes for 
which M(j) G ( 0 , 7 r) modulo 27r. 

A striking feature in the weak-coupling regime is the 
self-similar structure of the jc{4>) curve for ksT <C J: 
After the spontaneous M-fold increase of the unit cell, 
the chiral current exhibits a behavior that is similar to the 
one in the Meissner phase. For instance, for the VLi/ 2 j 
this is visible in Fig.|5][a) at a flux value of </> = tt and for 
VLs at p„ = 1/3 and 2/5 [see Fig. |5l[c) and (d)] at </> = 
27r/3 and 47 r/ 5 , respectively. As a result, the absolute 
values of the reversed current for (jj tt are several times 
larger than in the U = 0 case at the same value of cj) 
(compare the continuous black curve and the red dashed 
curve in Fig. |5l[a) near (p = ir). 

We studied the behavior of interacting bosonic par¬ 
ticles in the presence of a uniform magnetic flux on a 
two-leg ladder. Our work has three main results. First, 


we determined the range of stability of VLs predicted in 
00 and analyzed their microscopic structure. We ob¬ 
tained representative phase diagrams, showing that the 
system is very rich, also realizing the BLP phase 21\- 
Second, we observed the interaction-driven reversal of the 
chiral current in the vicinity of certain VLs. Third, we 
proposed an intuitive interpretation of the reversal via 
the spontaneous increase of the effective flux due to the 
enlarged unit cell in the VLs. We expect all these re¬ 
sults to influence experimental work on low-dimensional 
bosonic systems in the presence of gauge fields. 
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In this Supplemental Material we provide the definition of the chiral current, discuss its finite-size 
effects and we show additional data for its dependence on the interaction strength. We further 
discuss the phases mentioned in the main text in more detail, elaborate on the phase transitions 
and we show an additional phase diagram for density versus interaction strength. Additionally, we 
present patterns of the local current and density used for the characterization of the ground-state 
phases. Finally, we describe in more detail the methods used to obtain the results of the main text. 


SI. CHIRAL CURRENT 


The local current operators, defined through the stan¬ 
dard lattice version of the continuity equation (setting 
the lattice constant a = 1 and h = 1) for particle den¬ 
sities, are denoted by j| ^ (along the links of the ladder 
legs) and (on the rungs). Their explicit expressions 
are 

j^ r ~ 

jr = ,.a2,r - ■ (SI) 


In the thermodynamic limit {L —>■ 00 ), where contribu¬ 
tions from the rungs at the boundaries (r = 1 and r = L) 
of the ladder are negligible, or for periodic boundary con¬ 
ditions, the chiral current becomes 



r 



dEo 
d(j) ’ 


(S2) 


where Eq is the ground-state energy per particle [I| 
(in Secs. S1-S4 we will discuss ground-state properties 
(hence T = 0)). 

One can easily show that for non-interacting bosons 
{U = 0) and for < 2J sin | tan 


Jc — 


J? sin (j) 


8Jsin"' 2 V 1 sin 


(S3) 


For U = 0, the chiral current saturates at jc = Jsin^ 
in the Meissner-superfluid (MSF) phase, i.e., for J± > 
2J sin I tan 

In the complementary regime of strongly interacting 
bosons U ^ J, we also observe a quadratic dependence 
of jc on J± for small values of J±/J. In fact, such a de- 

S ndence holds for arbitrary values oi U/J for J± Jcj) 
. For strong hopping on the rungs J± J and for 
very strong interactions U/J ^ 00 , the chiral current 


can decay to zero in the Mott-insulating states even at 
half-integer filling. For instance, for p = 0.5, the chiral 
current is jc = sin^/(2J_L) [l|. Both in the U = 0 

and the G/J —>■ 00 regime, the particles circulate coun¬ 
terclockwise along the boundary given that 0 < cj) < n. 


S2. CHIRAL-CURRENT REVERSAL AND 
FINITE-SIZE DEPENDENCE OF THE CHIRAL 
CURRENT 


In Fig. 3 of the main text, we presented results for the 
chiral current for the largest system size (L = 160) used 
in our DMRG simulations. Strong finite-size effects are 
present in the vortex-liquid phase (VSF). In Fig. EH we 
further illustrate these finite-size effects and the current 
reversal effect for different fillings by showing additional 
results for the chiral current for systems with L = 80, 120 
and 160 rungs. By contrast, in the vortex-lattice phases, 
finite-size effects are practically invisible on the scale of 
the figure. 

As one can see in Fig. ISll upon increasing system size, 
the width of the small quasi-plateaux of almost constant 
chiral current present in the VSF phases diminishes, in¬ 
dicating that in the thermodynamic limit, the chiral cur¬ 
rent will become a continuous function of the flux in all 
phases as well as at the phase transitions between the 
vortex-liquid and the vortex-lattice phases presented in 
Fig. 3 of the main text. For a filling of p = 1 (presented in 
Fig. eh (d)), the transition from the vortex-lattice state 
with increasing U/J happens inside the Mott-insulating 
(MI) phase, with finite-size effects being much less severe. 

In the normalization adopted in Fig. ISll the value of 
the chiral current is independent of p and positive at 
17 = 0, implying that in all cases, the chiral current un¬ 
dergoes a sign change in the vortex-liquid phases as the 
transitions to the vortex-lattice states are approached. 
The data shown in the figure further shows that the 
chiral current reversal occurs irrespective of density or 
whether the underlying state is superfiuid (SF) or Mott- 
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FIG. SI. Finite-size dependence of the chiral current and dif¬ 
ferent fillings. We plot jc as a function of interaction strength 
U/J for the parameters of Fig. 3 of the main text (0 = O.Qtt, 
J± = 1.6J). The figure shows results for three system sizes 
[L = 80 (blue symbols), L = 120 (red symbols) and L = 160 
(black symbols)] and for particle densities (a) p — 0.4, (b) 
p = 0.5, (c) p = 0.8 and (d) p = 1. Due to the incommen¬ 
surate vortex density a significant system-size dependence is 
present in vortex-liquid superfiuid states (VSF), as opposed to 
the behavior in vortex-lattice states. For p = 1 we just show 
points for L = 120, which agree with the data for L = 80 on 
the scale of the figure. 


insulating. 


S3. GROUND-STATE PHASES 

Next we discuss details of the phase diagrams in the 
parameter plane for U/J = 2 and a non-integer 

filling as shown in Fig. 4 of the main text. Additionally, 
we present the phase diagram in the (p, C//J) plane for 
(j) = O.Ott and J± = 1.6 J in Fig. IS21 

It is quite remarkable that the ground-state phase di¬ 
agram shown in Fig. 4 of the main text is so rich, with a 
plethora of superfluid phases. Besides the MSF, the VSF, 
and the vortex lattice SF with vortex density pv = 1/2 
(VLi/ 2 -SF), a VL 1 / 3 -SF can get stabilized with a sponta¬ 
neously tripled unit cell. However, the stability domain 
of the VL 1 / 3 -SF [for the parameters of Fig. 4 of the main 
text] corresponds to flux values (f> > 27r/3 and there¬ 
fore, the increased (tripled) effective flux in the VL 1 / 3 -SF 
state does not result in a current reversal, confirmed by 
our simulations. For larger values of J±/J, an additional 
phase, predicted in mean-field approximation Q (see also 
additional recent work [^), with a spontaneously emerg- 



FIG. S2. Ground-state phase diagrams. Phase diagram in 
the (p, U) plane for (j) = O.Qvr and J± = 1.6J. The horizontal 
bold lines indicate Mott-insulating phases at commensurate 
fillings of p = 0.5,1, 1.5. 


ing imbalance in the particle numbers on the two legs 
exists, dubbed biased ladder phase (BLP). This phase, 
however, preserves the original unit cell (see Fig. lS5|) and 
thus does not trigger a chiral-current reversal either. 

In Fig. |S2l we present the phase diagram as a function 
of p and U/J for (/> = O.Ott, J± = 1.6J, p < 2 and U > 
J. The vortex lattice with p„ = 1/2 exists in a wide 
parameter range at sufficiently large densities (the critical 
lower density is H-dependent) and, for larger U/J > 3, 
is interlaced multiple times by the VSF phase as density 
increases. Moreover, there are Mott-insulating phases at 
densities p = 1/2,1, 3/2 [l|,[3-@|- The transition from the 
VL 1 / 2 -SF into the VSF phase induced by changing U or 
p is caused by a proliferation of soliton-anti-soliton pairs 
in the VSF state (see Fig.|S3]in the next section). 

The VSF is the only two-component Luttinger liq¬ 
uid [I| (central charge c = 2 ), while all other super- 
fluid phases are one-component Luttinger liquids (c = 
1). The vortex-lattice and Meissner phases are identi¬ 
fied by their characteristic local current configurations. 
They may be clearly discriminated from the vortex- 
liquid phase by calculating the central charge c extracted 
from entanglement-scaling properties [3,0 , similar to the 
strong-coupling U/ J ^ 00 regime previously studied in 
0 - 

All phase transitions in Fig. 4 of the main text and 
Fig.|S2]are second-order commensurate-incommensurate 
transitions except for those into the BLP phase, 
which are identified by a sharply increasing par¬ 
ticle density imbalance between the legs An = 
Sr (”1.^ “ ’^2,r) /N. 

There are also isolated Berezinskii-Kosterlitz-Thouless 
(BKT) da, [m critical points from Mott-insulating 
phases to superfluid states, along the paths of a constant 
particle density. Similarly the vortex lattice to vortex- 
liquid transitions for a constant vortex density are ex¬ 
pected to be of BKT type. 

The Mott-insulating region at p = 1 [shown also in 
Fig. ISlT dl] is richer and apart from a narrow Ulji/ 2 -mi 
region, a Meissner Mott-insulator state (M-MI) state is 
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realized with increasing U/J. Our simulations on finite- 
size ladders suggest the possibility of an intermediate 
vortex Mott-insulator (V-MI) region between the VL 1 / 2 - 
MI and M-MI phases. It remains yet to be confirmed 
whether this intermediate V-MI survives in the thermo¬ 
dynamic limit or whether a direct transition between the 
VL 1 / 2 -MI and M-MI phases is realized, similar to the 
case dX (j) = -K Bill- 


S4. LOCAL CURRENTS, PARTICLE 
DENSITIES AND CHIRAL CURRENT IN 
DIFFERENT PHASES 

Currents and densities atU = 2J, p = 0.8. The differ¬ 
ent superfluid phases discussed in the main text can be 
distinguished from each other by the structure of local 
particle currents and local particle densities. In Fig. [S3l 
we present typical examples for a cut through the phase 
diagram of Fig. 4 of the main text at J_l = 1.6 J. We show 
results for a central portion of the ladder with L = 120 
rungs. 

Upon increasing (j) (from bottom to top in Fig. [S3]), 
the system passes from the Meissner phase [MSF, 
Figs. [S3l' ell-fe3l] (stable for (jj < O.Ttt) to the vortex- 
lattice state VL 1 / 2 -SF (shown in Figs. |S3Kal)-(a3), real¬ 
ized for tt) by going through additional intermedi¬ 
ate vortex-liquid and vortex-lattice phases. As already 
mentioned, in the Meissner phase, currents are non¬ 
zero only along the ladder’s legs (and, for open bound¬ 
ary conditions, also on the boundary rungs r = 1 and 
r = L, not shown), hence the vortex density pv van¬ 
ishes, while the particle density is uniform in the bulk of 
the ladder. Adjacent to the Meissner phase, there is a 
vortex-liquid superfluid (VSF) with a finite vortex den¬ 
sity 0 < pv < 1/3 [shown in Figs. [S3Kdl)-(d3)]. For val¬ 
ues of O.TStt < < 0.827r, the vortex-lattice phase with 

Pv = 1/3 (VL 1 / 3 -SF) becomes the ground state [shown 
in Figs. [S3][cl)-(c3)]. Both local currents along the legs 
and along the rungs have a periodicity with a period of 
three times the lattice spacing, which is also exhibited by 
the local particle density. 

Upon further increasing the flux but still before the 
vortex lattice at p^ = 1/2 (VL 1 / 2 -SF) is stabilized, the 
system reenters into the vortex-liquid phase with vortex 
densities 1/3 < py < 1/2 [shown in Figs. [S3][bl)-(b4)]. 
In the vortex-liquid state and in the proximity of the 
VL 3 / 2 "SF, domain walls (or solitons) between the two de¬ 
generate ground states of the VL 1 / 2 -SF state, which we 
denote by |1) and |2), are clearly visible [see Fig. [S3][bl) 
and (b2)]. We call a soliton a defect in the periodic struc¬ 
ture of the VLi/ 2 "SF state with a |l)-like pattern to the 
left of the defect and a | 2 )-like pattern to the right side 
and anti-solitons a defect with a | 2 )-type pattern on the 
left side and a |l)-type pattern on the right side. The 
number of solitons and anti-solitons in the vortex-liquid 
state in the vicinity of the VL 3 / 2 "SF are identical, due to 
the topological confinement stemming from the double 


degeneracy of the ground state of the VL 1 / 2 -SF phase. 
The local density shows modulations near the positions 
of solitons and anti-solitons. In the VLi/ 2 "SF state, on 
the contrary, the local densities are constant in the bulk 
along the ladder and the local currents show a modu¬ 
lation with a period of two lattice cells. In general, in 
vortex-liquid phases, the particle density and local cur¬ 
rents exhibit a complicated structure, unlike the highly 
regular patterns of the vortex-lattice states. 

The maximal possible vortex density in this system 
can in fact not exceed py = 1/2. This limit is reached 
for values of the flux ^ tt, in agreement with the 
picture originating from the weak-coupling regime p 3 > 
1 , 0 [ll. The current configurations presented in our 
VL 1 / 3 -SF and VL 1 / 2 -SF phases, here obtained from our 
numerical simulations and shown in Figs. lS3( clVlc3l and 
Figs. [S3l' alV-('a3l. respectively, are different from those 
depicted in P,|^. Namely, and most importantly, around 
the vortices, particles circulate in the direction opposite 
to the behavior in the Meissner phase. Hence, the vor¬ 
tices that we observe resemble the vortices of type-II su¬ 
perconductors. We further note that in VL 1 / 3 -SF vortex 
lattices, the particle density is reduced at the corners 
of those plaquettes at which the vortices sit, which is 
not captured in the standard linearized bosonization ap¬ 
proach used in [l|,|3)0B- By contrast, in the VL 1 / 2 -SF, 
vortices occupy every other plaquette and hence there is 
no corresponding effect on the local particle density. The 
density pattern of the VL 3 / 3 -SF state is, however, in¬ 
homogeneous, yet nevertheless follows the periodicity of 
the local current configuration of the underlying vortex- 
lattice state: Vortices sit in regions of lower particle den¬ 
sity. In the VL 3 / 3 -SF, a periodically modulated local 
current and a periodic density modulation coexists with 
superfluidity. 

Unlike linearized bosonization, a mean-field ap¬ 
proach 0 yields a modulation of the particle density 
in vortex phases. However, in mean-field theory and 
at (^ = TT, the particle density is modulated from rung 
to rung with a period of two lattice constants, whereas 
our numerical simulations for fluxes close to tt (including 
4> = tt) produce VL 1 / 2 -SF states with a uniform particle 
density along the ladder. 

Chiral eurrent as a function of flux. The dependence 
of the chiral current as a function of flux <f) and as the 
system goes through the different ground-state configu¬ 
rations shown in Fig. [S3] is presented in Fig. IS4I For 
J± = I. 6 J, the chiral current does not undergo a sign re¬ 
versal in the VL 3 / 3 -SF state, since this state is stabilized 
only for (j) > 2 ttI3. A reversal occurs only near (j) it. 
Even though generically, the transition from Meissner 
to vortex-liquid phases is a second-order commensurate- 
incommensurate transition, in the numerical data for a 
finite system size shown in Fig. IS41 there is a finite jnmp 
in the chiral cnrrent. This jump, observed in our finite- 
size simulation, can be attribnted to the proximity of 
the multicritical point at which the Meissner superflnid 
phase, the vortex snperfluid and the BLP phase meet (see 
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FIG. S3. Configurations of local currents and densities in various phases. DMRG results for the cut at = 1.6J through 
the phase diagram presented in Fig. 4 of the main textt for O.Stt < < -k {U = 2J and p = 0.8). The configurations presented 

here correspond to the following flux values in descending order: (al) O.OStt, (al) O.OStt, (a3) O.Ott, (bl) O.SStt, (b2) O.SGtt, (b3) 
0.827r, (b4) O.TOtt, (cl) O.TTvr, (c2) O.TStt, (c3) 0 . 737 r, (dl) O.Tlvr, (d2) 0.697r, (d3) 0 . 687 r, (el) O.G4n, (e2) O.Gtt, and (e3) O.SGtt 
(L = 120 rungs). We only depict the behavior for the middle rungs of the ladder (40 < r < 80) in the bulk of the system. 
The local currents and local particle densities presented in Fig. 2(c)-(d) of the main text in the VL 1 / 3 -SF and VLi/ 2 "SF are 
cut from the central region of configurations (c2) and (a3) presented here. The length of the arrows encodes the absolute value 
of the local currents and the size of the circles encodes the onsite density. In the top four configurations (al)-(bl), the chiral 
current is reversed. 


Fig. 4 of the main text). 

Local currents and particle densities in the BLP phase. 
For completeness we also present results for the typi¬ 
cal behavior of local currents and particle densities in 
the BLP phase, obtained from DMRG simulations for 
L = 15, shown in Fig. [S5] (we choose a small system to 
illustrate the behavior on the boundary rungs r = 1 and 
r = L as well). The local currents are non-zero along the 
legs and edge rungs of the ladder and particles flow in 
counterclockwise direction. The absolute values of local 


currents (which for a finite system depend on the dis¬ 
tance from the edge rungs) are identical in both legs, 
while the density is imbalanced between the two legs. 
Thus, one can conclude that the behavior of the chiral 
current in the BLP superfluid is quite similar to that in 
the Meissner superfluid, where particles circulate along 
the ladder’s boundary, vanishing quickly on the rungs in 
the bulk since p„ = 0). Nevertheless, the BLP phase is 
distinct from the MSF since in the former there is a spon¬ 
taneously broken discrete Z 2 symmetry corresponding to 
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FIG. S4. Flux dependence of the chiral current. We display 
DMRG results for the chiral current along the cut at J± = 
1.6J through Fig. 4 of the main text {U — 2J and p — 0.8) 
for 0 < cf) < n, obtained for L = 120 rungs. 



FIG. S5. Typical pattern of local currents and local densities 
in the biased ladder phase (BLP). In the particular realization 
shown here the lower leg hosts the majority of particles. These 
DMRG results are for the parameters of Fig. 4 of the main 
text (p = 0.8 and U = 2J) and for a point inside the BLP 
phase {(j) = 0.87r, J±_ = 2.7J, L = 15). The length of the 
arrows encodes the strength of local currents and the size of 
the circles is proportional to the density. 


time reversal combined with exchange of the legs. 


S5. TEMPERATURE DEPENDENCE IN THE 
WEAK-COUPLING LIMIT 

In Fig. [S6l we present additional results on the tem¬ 
perature dependence of the chiral-current reversal. As 
mentioned in the main text one can observe the current 
reversal related to the presence of a VL 1/2 state close to 
()) = tt for ksT < J/2. The inset of Fig. [S6l shows the 
())-dependence of jc for several temperatures for Jj^ = J. 
For the other vortex-lattice phases the current-reversal 
effect is washed out by increasing temperature faster. 


S6. METHODS 

DMRG. For interaction strengths oi U/ J > 1 and par¬ 
ticle filling p < 2, we perform large scale numerical den¬ 
sity matrix renormalization group [T^ . [l8l| simulations 
which allow us to calculate zero-temperature properties 
for open boundary conditions and systems with up to 
L = 160 rungs, using 1000 DMRG states. Since the lo¬ 
cal Hilbert space of bosons is generally unconstrained we 
employ a cutoff of the maximal occupation. We typically 
use a cutoff of four bosons per site which is justified due 
to the repulsive nature of the onsite interactions and for 
U > J. By comparison with larger and smaller cutoffs 


we have verified the numerical accuracy of the quantities 
shown here. Figure 1571 shows the dependence of the chi¬ 
ral current on the bosonic cutoff rimax- For U > J, the 
results, obtained by keeping not more than rimax = 4 par¬ 
ticles per site, are already well converged. For 0 < U < J 
and by further increasing rimax, we have verified that the 
sign of the chiral current becomes positive again in the 
VSF phase for some U < Uc < J. In Fig. 3 of the main 
text we indicate this trend, based on data such as the ones 
shown in Fig. IS71 with a dashed line connecting DMRG 
data of ~ 1.5J and the known value at 17 = 0. More¬ 
over, the dashed lines in Fig. 3 of the main text indicate 
the kink that we expect to be present at the VL 1 / 2 -SF 
to VSF boundary, based on the data of Fig. [S7l 

Close to the VSF to VLp^-SF boundaries, the DMRG 
simulations tend to converge to metastable excited states 
with larger or smaller vortex density. Here, we have per¬ 
formed several calculations starting from different ran¬ 
domly chosen initial states. Selecting the lowest energy 
states gives the piecewise continuous chiral-current curve 
presented in Fig. [ST] and Fig. 3 of the main text. 

The BKT transitions between VL 1 / 2 -SF and VL 1 / 2 - 
MI at fillings p = 0.5 and p = 1.5 shown in Fig. IS2I and 
Fig. 3 (a) of the main text are determined by numeri¬ 
cally analyzing correlation functions. At the transition 
from the VL 1 / 2 -SF to the VL 1 / 2 -MI, bosonization pre¬ 
dicts that the single-particle correlation function decays 
as {a\.^ai^r+x) (up to logarithmic corrections). 

The transition point can be extracted accurately from 
the finite-size scaling behavior of peaks in the quasi¬ 
momentum distribution function (l9l| . 

For the BKT transition from the superfluid to the 
Mott insulator at unit filling p = 1, which happens in¬ 
side the vortex lattice at p„ = 1/2 state presented in 
Fig. [S2] and Fig.[Sl](d), it can be shown using bosoniza¬ 
tion that the single-particle correlation function decays as 



T/j 

FIG. S6. Temperature dependence of the chiral current in the 
weak-coupling limit. The coloring in the main panel encodes 
the strength of the chiral current jc/J as a function of tem¬ 
perature ksT/J and rung-coupling ,J±/J for (f = 0.927r. The 
curve defined as jc/J = 0 is highlighted with a dotted line. 
The inset shows jc/J for Jx/J = 1 as a function of the flux 
(j) for several temperatures ksT/J = 0.01 (dotted blue line), 
0.1 (dashed red line) and 0.2 (solid black line). 
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FIG. S7. Dependence of the chiral current on the bosonic cutoff Umax- DMRG calculation for the parameters of Fig. 3, (j) = O.Qtt 
and L = 40 rungs for density (a) p = 0.5, and (b) p = 0.8. 


{op j.ai,r+x) ^ x~» (up to logarithmic corrections). Other 
phase transitions are estimated from the behavior of the 
chiral current as well as from the local-current structure 
(such as the transitions indicated in Fig. [ST] and Figs. 3 
of the main text). 

Transfer-matrix approach. For large particle densities 
per ladder site p 1, we follow an approach used in the 
theory of Josephson-junction ladders, where each super¬ 
conducting island is proximity-coupled to its three neigh¬ 
bors. It is convenient to introduce a density phase rep¬ 
resentation of bosonic operators aj ^ = yJp-\- 
where 5ppr = — P and describe the density fluc¬ 

tuations and the phase of the condensate at site (Z,r), 
respectively. For U Jp, neglecting the charging effects 
of the superconducting islands [3) [3 j the Hamiltonian 
of Eq. (1) in the main text becomes 

L 

H —>■ —2Jp ^ cos(0^,r-|-i — Ol^r) 

t=l,2-,r=l 

L 

-2J^p^ cos(0i,r - 6*2,r -(S4) 

r—1 


with local variables 9pr S [0, 27r). In Eq. (IS4I1 . it is more 
convenient to impose periodic boundary conditions. The 
transfer matrix P, connected with the partition function 
Z = tr was obtained in [I^. Eollowing 0 , we nu¬ 
merically calculated the eigenvalues A„ of the transfer 
matrix (where n = 0,1,2... and |A„| > |A„+i|) in the 
thermodynamic limit. The periodicity of the argument of 
the second largest eigenvalue as a function of flux yields 
the vortex density Pvif’) = Arg[Ai(</>)]/27r for different 
temperatures 0 presented in Fig. 5 (b) in the main 
text. The chiral-current curves corresponding to differ¬ 
ent temperatures, presented in Fig. 5 (a) in the main 
text, are obtained from 


jcif) = - 


ksT 91n Z 
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